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I. INTRODUCTION 



Exclusive nonleptonic B meson decays provide information of important dynamics such as penguin contribu- 
tions, CP asymmetries, and weak and strong phases. However, it is difficult to analyze these processes because 
of their nonperturbative origin. To simplify the analysis, the factorization approximation (FA) has been applied 
to decay amplitudes, under which nonfactorizable and annihilation contributions are neglected, and factorizable 
contributions are expressed as products of Wilson coefficients, meson decay constants, and hadronic transition 
form factors. In this paper we shall compare three approaches to nonleptonic charmless B meson decays, which 
are more or less related to FA: the FA approach pj], the perturbative QCD (PQCD) approach and the 

Beneke-Buchalla-Neubert-Sachrajda (BBNS) approach ||||. We take the B — > nir decays as an example, and 
explain how factorizable, nonfactorizable, and annihilation contributions, infrared-cutoff and scale dependences, 
and strong phases are treated in these approaches. 

The effective Hamiltonian for the B — > itit decays is given by [Q 
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C 1 { l i)O i «\v) + C 2 {tiO%\») + Y,Ci{ii)O i { l x) , (1) 



v q—u,c L 3 

with the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements V q = V* d V g b and the four-fermion operators 
= {d t qj) V -A{q 3 bi) V -A , of = {diqi) V -A{qjbj)v-A , 



3 = (dibi) V -A^2{<ljQj)v-A , 4 = {dib ) v -A^iqjqdv -A , 

q q 

5 = (dibi)v-A*^2(qjqj)v+A , 6 = {djbj)y- A /MjQi)v+A , 



°t = 2(dibi)v-A^2e q (<ijQj)v+A , O s = - {d % b ) v _ A ^2e q (q J q i ) v +A , 
9 q 

Og = -(d l b i ) v ^ A ^2e q (q :i qj) V - A , O 10 = -(dibj) V -A ^2 e q%Qi)v-A , (2) 
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where i, j are color indices and /j, a renormalization scale. 
Decay amplitudes are written as 

M cx (7t(P 2 )7t(P 3 )\H cS \B(P 1 )) , 

oc ^C i (M)(7r(P 2 )7r(P 3 )|Oi(M)| J B(Pi)) , (3) 

i 

with the meson momenta Pj. For convenience, we assume Pi = (P^ - , Pj - , Pit) = (-^b/v / 2)(1 ) IjOt), 
P2 = (Mb/\/2)(1, 0, 0t), and P3 = (A/b/\/2)(0, 1, 0t), Mb being the B meson mass. The dependences on 
the renormalization scale /u, cancel between the Wilson coefficients Ci{n) and the hadronic matrix elements 
(7T7r|Oi(/i)|P), such that the decay amplitudes are scale-independent (and also scheme-independent). In the 
practical derivation of the effective Hamiltonian, an infrared cutoff for loop integrals needs to be introduced, 
since four-fermion amplitudes may contain soft divergences, which are related to nonperturbative dynamics 
involved in the matrix elements. If the infrared cutoff appears in the form that external fermions are off-shell, 



the derivation will even become gauge dependent gj. These potential infrared-cutoff and gauge dependences 
are usually hidden in the matrix elements. 

In the FA approach the matrix element of, say, O2 is expressed as 

( 7 r(P 2 )7r(P 3 )|02(M)|5(P 1 )) « (7r(P 2 )|(d igi )v_A|0)(7r(P 3 )|fe6,-)v-A|S(Pi)) , 

ocf w F B «(q 2 = MZ), (4) 

where / w is the pion decay constant and F B7T (M%) the B — > 7r transition form factor evaluated at g 2 = M 2 , M„. 
being the pion mass. The form factor itself is parametrized by models. Though analyses are simple under FA, a 
serious problem arises. When decay amplitudes are expressed as the products of /^-dependent Wilson coeficients 
and /i-independent hadronic form factors, they depend on the scale \x. A cure of this scale dependence has been 
proposed in ||[l0]], which will be discussed in Sec. II. 

The PQCD approach to exclusive B meson decays has been developed some time ago 11 1^] . PQCD is a 



method to separate hard components from a QCD process, which are treated by perturbation theory. Non- 
perturbative components are organized in the form of hadron wave functions, which can be extracted from 
experimental data. This formalism, so-called factorization theorem, has been successfully applied to various 
semileptonic decays 0,^2) and nonleptonic (charmed and charmless) decays [|[| 13 1^] . For the B — ► tttt 



decays, the hard amplitude involves six external on-shell quarks, four of which correspond to the four-fermion 
operators and two of which are the spectator quarks in the B and 7r mesons. Since nonperturbative dynamics 
has been factored out, one can evaluate all possible Feynman diagrams for the six-quark amplitude straight- 
forwardly, which include both factorizable and nonfactorizable contributions. Factorizable and nonfactorizable 
annihilation diagrams are also included. That is, FA is not necessary in the PQCD approach. Without FA, the 
scale independence can be achieved easily [fl7f . It has been shown that nonperturbative meson wave functions, 
being universal, are the same for various topologies of diagrams in all decay modes, which contain the B and 7r 
mesons. 

Recently, Beneke et al. proposed an alternative formalism for two-body charmless B meson decays [^) . They 
claimed that factorizable contributions (the form factor F B7V in this case) are not calculable in PQCD, but 
nonfactorzable contributions are in the heavy quark limit. Hence, the former are treated in the same way as in 
the FA approach, and expressed as products of Wilson coefficients and F Bt! . The latter, calculated as in the 
PQCD approach, are written as the convolutions of hard amplitudes with three (P, tt, 7r) meson wave functions. 
Annihilation diagrams are neglected as in FA. Therefore, this formalism can be regarded as a mixture of the 
FA and PQCD approaches. Values of form factors at maximal recoil q 2 — M 2 and nonperturbative meson wave 
functions are all treated as free parameters. The BBNS approach then involves more parameters compared to 
FA and PQCD. This approach has been applied to the B — > irir and Ktt decays ||. 

In this work we shall make a detailed comparision of the above three approaches to nonleponic charmless B 
meson decays. We discuss their treatments of the infrared-cutoff and scale dependences in Sec. II, factorizable 
contributions in Sec. Ill, nonfactorizable contributions in Sec. IV, and annihilation contributions and strong 
phases in Sec. V. It will be shown that the CP asymmetry of the B^ — > tt ± tt t decays predicted in PQCD is 
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much larger than in the FA and BBNS approaches. Hence, the measurement of CP asymmetries will provide an 
appropriate test of these predictions. Section VI is the summary. Note that final-state- interaction (FSI) effects 
are always assumed to be absent in the above three approaches. The neglect of FSI effects has been argued in 



II. SCALE DEPENDENCE 

As stated in the Introduction, a serious problem of FA is the dependence of physical predictions on the 
renormalization scale fi. A plausible solution to the aforementioned problem has been proposed in |) 10 1. The 



idea is to extract a /i-dependent evolution factor <?(/i) from the matrix clement (tttt\0(i-i)\B) , and combine it 
with the /i-dependent Wilson coefficient C(/i): 

M (x C(n){mr\0(ji)\B) = C{^)g{p){^\0\B) trcc = C eS (inr\0\B) tIee , (5) 

where the effective Wilson coefficient C is ^-independent. The extraction of g(/j.) involves the mixing of 
effective operators below the scale y, and above the lower scale mi,, mt, being the b quark mass. After making a 
physical amplitude explicitly scale- independent, FA is applied to the matrix element at the tree level, 

(7T7r|0|S) trcc oc UF B * . (6) 

This strategy has also been adopted in the BBNS approach. 

The above issue is in fact much subtler. One-loop corrections to on-shell b quark decay amplitudes are 
usually infrared divergent ||]. To remove these infrared divergences, the conventional treatment is to assume 
that external quarks are off-shell by — p 2 . This assumption, however, results in a gauge dependence of one- loop 
corrections. The evolution factor </(/i) then depends on an infrared cutoff —p 2 and on a gauge parameter, which 
are originally implicit in the matrix element. Therefore, the modified FA, though removes the scale dependence 
of a decay amplitude, often introduces the infrared-cutoff and gauge dependences. 

The controversy of the scale, infrared-cutoff and gauge dependences can be resolved in the PQCD approach 
[fl7f . In this formalism partons, i.e., external quarks, are assumed to be on shell, and both ultraviolet and infrared 
divergences in radiative corrections are isolated using the dimensional regularization. With external quarks being 
on shell, gauge invariance of a decay amplitude is guaranteed under radiative corrections. The ultraviolet poles 
are subtracted in a renormalization scheme. The infrared poles, which correspond to nonperturbative dynamics 
of decay processes, are absorbed into meson wave functions. How much finite piece of radiative corrections is 
absorbed into wave functions along with the infrared poles depends on a factorization scheme, or equivalently, 
on a factorization scale /i/. In the PQCD approach fj,f has been chosen as the inverse of the meson transverse 
extent, 1/6. In other words, dynamics below (above) the factorization scale is absorbed into a nonperturbative 
wave function (perturbative hard amplitude) , and a decay amplitude is factorized into convolution of these two 
subprocesses. 

For nonleptonic B meson decays, contributions above the factorization scale involve two scales: Myy, at which 
the matching conditions of the effective weak Hamiltonian to the full Hamiltonian are defined, and the typical 
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scale t of a hard amplitude, which reflects specific dynamics of a decay mode. The scale t, corresponding to the 
lower scale mj, in the modified FA, is chosen as the maximal virtuality of internal particles in a hard amplitude. 
The analysis of next-to-leading-order corrections to the pion form factor has indicated that this choice 
minimizes higher-order corrections. 

Without assuming FA, the fi dependence of a decay amplitude is calculated based on six-fermion hard ampli- 
tudes in the PQCD approach. Radiative corrections produce various types of logarithms: \n(Mw / 'fj,) , ln(/i/i) 
and ln(£6). The renormalization-group (RG) summation of the first type of logarithms leads to the Wilson 
coefficients C(y) in Eq. (||). The second and third types are summed into the evolution factor 

g(fi)= gi (ti,t)g 2 (t,l/b) . (7) 

The above expression implies that g in fact contains two pieces, g\ and g 2l which are governed by different 
anomalous dimensions. The anomalous dimension of g\ is the same as of C(/i). Below the hard scale t, loop 
corrections associated with spectator quarks contribute, such that the anomalous dimension of g 2 differs from 
that of .91 0. 

Insert Eq. (f7|) into Eq. (||), the effective Wilson coefficient is written as 

C cS = C(t)g 2 (t,l/b), (8) 

which consists of two stages of RG evolutions. We have extended the Wilson evolution directly from Mw down 
to t without further considering the mixing of effective operators. The scale independence of a decay amplitude 
is thus constructed explicitly. The Wilson evolution is then followed by the second-stage evolution from t down 
to 1/6, which plays the role of the infrared cutoff — p 2 . 

Together with the second-stage RG evolution, there also exists the Sudakov evolution exp[— s(P, b)] from the 
resummation of double logarithms In 2 (Pb) [jl9|,^o| , which arise from the overlap of collinear and soft divergences 
in radiative corrections to meson wave functions, P denoting the dominant light-cone component of meson 
momentum. The effect of the Sudakov evolution will be explained in the next section. After summing all 
various large logarithms, the hard amplitude H (t) can be evaluated at the characteristic scale t perturbatively. 
H(t) contains all possible factorizable and nonfactorizable diagrams. 

At last, we address the treatment of the dependence on the factorization scale in the PQCD approach. Both 
the hard amplitudes and the meson wave functions <p(x,b), x being the momentum fraction associated with 
a valence quark, depend on factorization schemes. For example, 4>{x, b) and 4>{x, 6/2) are different. However, 
a decay amplitude, as convolution of the two subprocesses, does not depend on factorization schemes. The 
wave functions, though not calculable, are universal (process- independent), since they absorb long-distance 
dynamics, which is insensitive to the specific decay of the 6 quark into light quarks with large energy release. 
In the practical application 4>{x, b/c) for the factorization scale c/6, c being a constant of order unity, can be 
determined from experimental data of some decay modes. These wave functions are then employed to make 
predictions for other modes, whose factorizations are constructed at the same scale c/6. With this prescription, 
the dependence on the factorization scale is removed. Usually, c is chosen as unity. 
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A three-scale factorization formula for exclusive nonleptonic B meson decays is then written as 

C(t) <g) H(t) <g) (/>(x, b) ® g 2 (t, b) <g) exp [s(P, b)} . (9) 

Note that Eq. (^) is a convolution relation, with internal parton kinematic variables x and b integrated out (t 
depends on x and b as defined later). All the convolution factors, except for the wave functions (f>(x,b), are 
calculable in perturbation theory. Comparing Eq. (||) with Eq. (^|) , the tree-level matrix element is read off as 

(7T7r|0|£) trcc = H(t) <g> 4>{x, b) ® exp [-s(P, 6)] . (10) 

If choosing t as mj, Eq. (||) reduces to a simple product of the constant Wilson coefficient C(m(,) and the 
hadronic matrix clement. 



III. FACTORIZABLE CONTRIBUTIONS 



In the FA approach factorizable contributions are expressed as products of Wilson coefficients and transition 
form factors, which are then parametrized by some models. This treatment of factorizable contributions has 
also been adopted in the BBNS approach. It has been argued that the B — > ir transition form factor F B ~* can 
not be evaluated in PQCD, because the on-shellness of internal particles in the hard amplitude is not suppressed 
by meson wave functions For example, the lowest-order hard amplitude in Fig. 1 consists of the internal b 
quark and gluon propagators proportional to 

I = I = 1 (11) 

[(Pi-fc 3 ) 2 -rn2](fc 1 -fc 3 )2 (2P 1 -k 3 )(2k 1 -k 3 ) x x x%M B ' V ' 

where k% (k 3 ) is the momentum of the spectator quark in the B (tt) meson, and x\ = k^ / {x 3 — k^ / P 3 ~) the 
corresponding momentum fraction as indicated in Fig. 1. To derive the above expression, we have neglected the 
mass difference A = Mb ~ mi,. Obviously, the pion wave function, if proportional to x 3 (l — x 3 ), does not remove 
the power divergence at x 3 — > from the on-shell internal particles. The assumption that internal particles in 
a hard amplitude should be off-shell is then violated, and PQCD is not applicable to F BlT . 

A different viewpoint has been adopted in the PQCD approach. Since the end-point divergence is not of 
the pinched type, which is absorbed into a wave function, we argue that it can be either killed by a pion wave 
function vanishing faster than 13 as 13 -> 0, or smeared out by parton transverse momenta kj- pl[ |. Including 
kr, Eq. (^lj) becomes 

(x 3 M B + k 2 3T )[x lX3 M B + (k 1T - k 3T ) 2 ] ' (12) 

which has no singularity at x 3 — * 0. As argued in ]2l| , partons can carry small amount of transverse momenta. 
When the end-point region with x — > is not important, it is proper to approximate Eq. ( |l2| ) for the hard 
amplitude by Eq. (|ll]), and the transverse degrees of freedom of partons in meson wave functions can be 
integrated out. When the end-point region is important as in the case of F B7V , one should keep kr for a 
consistent analysis. Once the transverse degrees of freedom of partons are taken into account, the factorization 
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of a QCD process should be performed in the 6 space, 6 being the variable conjugate to far |P0| . This is how 
the factorization scale 1/6 is introduced in the previous seciton. 

As stated before, the scale t should be chosen as the maximal virtuality of internal particles in a hard 
amplitude, 



according to Eq. (|lj). The Sudakov factor exp[— s(P, b)] suppresses the long-distance contributions in the large 
b region, and vanishes as b = 1/Aqcd, Aqcd being the QCD scale. This factor then guarantees that most 
contributions to decay amplitudes come from the region with small coupling constant a s (t) It has been 

shown that 97% of the contribution to F Bn arises from the region with a s {t)/Ti < 0.3 [Q. It indicates that F Bw 
at large recoil of the pion is well within the perturbative regime, and that dynamics from hard gluon exchanges 
dominates. This is the reason PQCD with Sudakov suppression is applicable to exclusive decays around the 
energy scale of the B meson mass The above conclusion is consistent with the following physical picture. 
The B — > 7r transition occurs via the b — > u decay with the outgoing u quark being energetic. The spectator 
quark of the B meson is more or less at rest. The probability that this spectator quark and the u quark with 
large relative velocity form the pion is strongly suppressed by the pion wave function and by the Sudakov factor. 
Therefore, the process dominates, where a hard gluon is exchanged so that the two quarks are aligned to form 
the pion. 

We need to illuminate another concept about the application of PQCD. The smallness of PQCD predictions 
compared to experimental data does not definitely imply that perturbative contribution is not important p3]. 
PQCD predictions depend on wave funcitons one adopts. For example, a smooth B meson wave function 
decreases PQCD predictions for F Bt! by a huge extent. However, this smooth B meson wave function may 
not be correct, such that a conclusion drawn from it is not valid. In our analysis the B meson and pion wave 
functions are determined from the data of the pion form factor and of the B — > Dir(p) decays. We then employ 
the extracted wave funcitons to predict F B7r and examine the dominance of perturbative contributions. In this 
way the ambiguity in choosing wave functions and the model dependence of PQCD predictions are reduced. 

Below we highlight the enhancement of penguin contributions observed in the PQCD approach, and its role 
in the explanation of the B — > rnr and B — > Kir data. For simplicity, we present the observation by means of 
the FA approach. Consider the ratios R and R^, 



t = max (^1X3 M B , 1 /i 3 "M B , 1/6) , 



(13) 



R 



Br( J B° -> AT ± 7r = F) 



a\ + 2a K X 2 R b cos </) 3 



(14) 



R 



IT 



Br(B± -> AT°7rF) ' 
Bi(B° -> A^ti-t) 
1 Br(_B° -> t^tt^) 



\ 2 R b [R b + 2a^{R b - cos0 3 )] 



a 2 K + 2dK^ 2 R b cos (f) 3 



(15) 



with the Wolfenstein parameters A and R b , and the Wilson coefficients 



ar = C 2 + -, 
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N c 2 q V iV, 
N c being the number of colors. The factor 

a 4 + 2r KM a 6 _ M k(tt) 

y ' ai y > M B {m s(u) +m d ) 

being negative, represents the ratio of the penguin contribution to the tree contribution in the Ktt^tt) mode, 
where M K ^) is the kaon (pion) mass and m q , q = s, u, and d, the q quark masses. 



It is obvious that the data R ~ 1 |23| imply the unitarity angle 4>3 ~ 90°, no matter what ax, A and Rb 
are. Therefore, R is an appropriate quantity for the determination of 03. While to determine <f>^ from the data 
of the ratio R^ ~ 4 [ p3| , one must have precise information of o,k , o>-ki A and Rb- It can be shown that the 
extraction of 4>3 from R^ depends on these parameters sensitively. Hence, R^ is not an appropriate quantity 
for the determination of fa. To explain R^ ~ 4, a large |a ff | ~ 0.1 corresponding to = 2m„ = 3 MeV and a 
large fa ~ 130° must be adopted in the FA approach [§3. This is obvious from Eq. (|l5|), since a large ^3 leads 
to a constructive interference between the first and second terms in the numerator of R n . In the modified FA 
approach with effective number of colors N^ s , a large unitarity angle ^3 ~ 105° is also concluded |2fJ. 

It is interesting to ask whether one can explain R v ~ 4 using reasonable quark masses m c i = 2m u ~ 10 
MeV and a smaller ^3 ~ 90°, if the extraction from R is given a higher weight. The answer is positive in 
the PQCD approach. An essential difference between the FA and BBNS approaches and the PQCD approach 
is that transition form factors can be calculated in the latter as argued above. Therefore, we do not assume 
that the form factors associated with the Wilson coefficients 01,4 and with are the same. An explicit PQCD 
calculation indicates that the ratios clk.-k are more complicate, which reduce to Eq. ( p"7| ) only in some extreme 
kinematic conditions of partons B. Moreover, the choice of the hard scale t, at which Wilson coefficients are 
evaluated, is also different. It is chosen arbitrarily as mj or m&/2 in the FA and BBNS approaches, but as 
virtuality of internal particles of a hard amplitude in PQCD, such that evolution effects vary among the above 
form factors. 

It has been noticed that <xi(/x), a 4 (^) and ae(p) have dramatically different /i dependences for /1 < mf,/2: |ct4| 
and \ae\ increase with the decrease of /j, faster than a\. The choice of the scale in PQCD makes possible that t 
runs to below rrib/2, and that penguin contributions are enhanced over tree contributions. It has been observed 
that the evolution effect enhances the ratios |aK,7r| by about 50% |Q. Hence, even adopting = 2m u ~ 10 
MeV, we can account for R n ~ 4 using ^3 ~ 90°. That is, the data of R n do not demand large <j>^. We 
emphasize that the hard scale t with an average 1.4 GeV, corresponding to a s (t)/n ~ 0.13, is still large enough 
to guarantee the applicability of PQCD [Q. 

IV. NONFACTORIZABLE CONTRIBUTIONS 

Nonfactorizable contributions arc neglected in the FA approach. However, it has been known that this 
approximation is not always appropriate, especially for the decay modes, whose factorizable amplitudes are 
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proportional to the Wilson coefficient a 2 = C\ + C 2 /N c . Since a 2 is small, factorizable contributions and 
nonfactorizable contributions are of the same order, and the latter become important. This is the reason it is 
difficult to resolve the controversy in the branching ratios of the B — > J/ipK^ decays in FA The conventional 
way to introduce nonfactorizable contributions is to vary number of colors in the Wilson coeficients a^. The 
resultant effective number of colors, N% , being process-dependent parameters, render the FA approach less 
predictive. 

To go beyond FA, it has been proposed to evaluate nonfactorizable amplitudes for two-body charmless B 
meson decays using PQCD factorization theorem in the BBNS approach ||. The motivation is that the power 
divergence appearing in F B7r does not exist in nonfactorizable amplitudes, since soft divergences cancel between 
the two diagrams in Fig. 2 J5J. Hence, these diagrams contribute to the hard parts of the B — > irir decays, 
and are calculable in PQCD. We emphasize that the above soft cancellation is not sufficient to justify the 
applicability of PQCD. No matter what value Mb is, one can always obtain the soft cancellation between the 
two nonfactorizable diagrams. However, it is obvious that PQCD is not applicable to B meson decays if Mb 
is as low as 1 GeV. The criterion for the applicability of PQCD is stronger: PQCD is applicable only when 
higher-order corrections to decay amplitudes are under control. 

We shall show explicitly, following the reasoning in ]2^| , that the simple factorization formula in Jj| for a 
nonfactorizable amplitude suffers large higher-order corrections. The expression is 

(18) 

where the constant C is related to a color factor, the denominator xix comes from the hard gluon propagator, 
and another denominator x comes from the internal quark propagator. If the hard part is characterized by Mb, 
the coupling constant a s evaluated at the scale fi ~ Mb is small and higher-order corrections are negligible. 
However, it is not the case. When including higher-order corrections that cause the running of a s , the logarithm 
a s (n) \n.{x\xM B / /i 2 ) occurs, with X\xM B being the invariant mass of the hard gluon. In order to diminish this 
logarithmic correction, the appropriate choice of the scale is /i — ^/x\xMb- The coupling constant a s in Eq. (|l8| ) 
becomes running and should move into the convolution relation, leading to 

/dx dx dx 

— — -4>tt{x 2 ) — -(j)„(xs)a s (-Jxix 3 M B ) ■ (19) 
Xl x 2 x 3 

Assume the following models for the meson wave functions: 

(f> B (x) 

tf B {x) 

The B meson wave functions 4>b 0J and 4>' B possess peaks at small x that depend on the shape parameters 
lob and to' B , respectively. According to j5|, the values of the shape parameters can be fixed via the parameter 
X B = 0.3 GeV: 



/ = Ca s [ — <Pb(xi) [— faix) 

J Xl J X 



N B \/x(l - x) exp ^~ 
N' B x 2 (l -x) 2 exp 
-^-f v x(l - x) . 



M%x 2 
M 2 B x 2 



(20) 
(21) 
(22) 
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J ° i J ' , = -j-^ = 17-6 , (23) 

which gives lub — 0.65 GeV and uj' b = 0.25 GeV. The normalization constant Nb is related to the B meson 
decay constant fs through L (f>B{x)dx = fg / (2^/2N c ). The asymptotic pion wave function in Eq. (22) has 
been adopted in J5j. 

Using the wave function 4>b{4>'b)i we nn d t na ^ only about 20% (30%) of the full contribution to the nonfactor- 
izable amplitude comes from the perturbative region with a s /7r < 0.3. To derive the above percentages, we have 
chosen Aqcd = 250 MeV for the running a s , and frozen the scale Jx{ x 3 Mb to 0.4 GeV, as it goes to below 0.4 
GeV. We conclude that Eq. (|l9| ) is dominated by soft contribution from the end-point region with x%,x — > 0, 
i.e., with a s > 1, where the perturbative expansion is not justified. Therefore, nonfactorizable contributions 
can not be evaluated reliably in the BBNS approach. This conclusion is consistent with that drawn in J2(J: 
PQCD without Sudakov suppression is not applicable to exclusive processes for an energy scale below 10 GeV. 
The formalism in || is self-consistent as Mb — > oo, but not in the case with Mb ~ 5 GeV. 

In the end-point region where soft contributions are important, nonfactorizable amplitudes should be treated 
in a more careful way. That is, the dependence on the parton transverse momentum fc^ is essential, and should 
be included as argued in the previous section. The argument of a s is also set to t, the maximal virtuality of 
internal particles of a hard amplitude, which is similar to that in Eq. (|l9|), but depends on the additional scale 
1/6. a s (t) is then always small even in the end-point region, since 1/6 is large under Sudakov suppression |2lf| . 
It has been shown that almost 100% of full contributions to nonfactorizable amplitudes arise from the region 
with a s (t)/n < 0.3. That is, the evaluation of nonfactorizable contributions is reliable in the PQCD approach. 

There is another difference between the BBNS and PQCD approaches in the treatment of nonfactorizable 
contributions. The momentum of the light spectator quark in the B meson has been ignored in the former jjj], 
such that quark propagators in the hard part always remain time-like: 

1 1 _ 1 

(k 2 + k 3 - fcx) 2 + ie ~ (k 2 + fc 3 ) 2 ~ x 2 x 3 Ml ■ ( } 

This approximation drops the strong phase from the kinematic region where the virtual quark becomes on-shcll, 
(k 2 + k 3 — ki) 2 ~ 0. In the PQCD approach nonfactorizable contributions are complex: 

1 



{k 2 + k 3 ~ h) 2 + ie 



P 



tS [{k 2 + k 3 - h) 2 ] , (25) 



_(k 2 +k 3 - ki) 2 _ 

where P denotes the principle- value prescription. As shown in fllq] , the above imaginary part is not the main 
source of strong phases for the B — > tttt decays, because nonfactorizable contributions are only few percents 
of factorizable ones. However, it will become important in decay modes whose factorizable amplitudes are 
proportional to the small Wilson coefficient a 2 , or absent, such as the B® — > K^K T decays pS{] . 

V. STRONG PHASES 

In the FA and BBNS approaches annihilation contributions have been neglected. Though annihilation am- 
plitudes from the operator O4 for the B — > nn decays indeed vanish because of helicity suppression, those from 
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Oq do not. As shown in Eqs. (Bl) and (B6) of j29f| . there is a large enhancing factor 

Ml 



(m b + m d ){m u + m d ) 

associated with annihilation contributions from Oq, which is much larger than the one 

Ml 



(26) 



(27) 



(m b + m d ){m u + m d ) 

associated with usual penguin contributions. Though annihilation contributions are suppressed by a factor 
/b/Mb pej] , they can be comparible with penguin contributions. It has been observed that the annihilation 
diagrams Figs. 3(e) and 3(f) contribute large imaginary parts (strong phases) in the PQCD approach 0]. 

In the approaches where annihilation contributions are neglected, strong phases must be introduced via 
the Bander-Silverman-Soni (BSS) mechanism shown in Fig. 4(a) However, this mechanism, contributing 

through a charm quark loop, is suppressed by the charm mass threshold f§: 

q 2 = (x 2 P 2 + x 3 P 3 ) 2 = x 2 x 3 M% > Am 2 c , (28) 

m c being the c quark mass. Obviously, the large x 2 and x 3 regions are not favored by the pion wave funtion 
in Eq. (^2|). . Moreover, the diagram with a charm quark loop is of higher order in the PQCD formalism. An 
exact numerical analysis has indicated that the imaginary contribution from the BSS mechanism is smaller than 
that from the annihilation diagrams by a factor 10 Q. 

Another source of strong phases in the FA and BBNS approaches is the diagrams in Fig. 4(b), whose imaginary 
contributions appear in the effective Wilson coefficients in Eqs.(4)-(8) of H], for instance, 

n u - r + ° 2 4- " s Cf r f 

F = -121n-^-18 + # + # J , (29) 

with 

Si = -0.5 - 3ttz , (30) 



and fl 1 being real. We emphasize that Eq. (30) is derived from a configuration, in which the outgoing u quark 
carries the full pion momentum. Again, this configuration is strongly suppressed by the pion wave function and 
by the Sudakov factor. That is, the imaginary contribution from the above source has been overestimated. 

Below we make a quantative comparision among the FA, BBNS and PQCD approaches. Figure 3 collects 
all the Feynman diagrams for the hard amplitude of the B d — > tt + tt~ decay, including factorizable and non- 
factorizable tree, penguin and annihilation diagrams, whose contributions are listed in Table I. The results 
obtained in the BBNS approach are also shown. The FA results are the same as the BBNS ones, but with the 
nonfactorizablc amplitudes neglected. We observe 

1. The tree contributions are almost equal in the three approaches, but the penguin contributions are larger 
in PQCD, corresponding to a larger ratio \a^ \ defined by Eq. ( p^ ) (about twice of those in FA and BBNS), 
i.e., corresponding to the dynamical penguin enhancement. 
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2. The tree and penguin contributions are real in PQCD, but complex in FA and BBNS due to the complex 
effective Wilson coefficients, which arc the main source of strong phases in FA and BBNS. 

3. The annihilation contributions are neglected in FA and BBNS, but of the same order as the penguin 
contributions in PQCD, which are the main source of strong phases in PQCD. 

4. Nonfactorizable contributions are complex in PQCD but real in BBNS, as explained in the previous 
section, and neglected in FA. Note that different B and it meson wave functions have been employed as 
deriving the results in Table I [fi]J|] . Nonfactorizable contributions are destructive relative to the penguin 
contributions, if concering only the real parts, and negligible for the B — » ttt: branching ratios in both 
approaches. 

5. The dominant nonfactorizable contribution M p in BBNS is due to the sharp hard parts in the end-point 
regions of momentum fractions x as stated in Sec. III. 

6. The total tree and penguin contributions, including both factorizable, nonfactorizable and annihilation 
ones, are 

Tpqcd = (71.87 + 3.2(H) x 1(T 3 , P P q C d = (-7.47 + 4.64z) x 1(T 3 , 

Tbbns = (74.93 + 1.13*) x 1CT 3 , P BB ns = (1.13 - 1.27*) X 1(T 3 , 

T FA = (75.70 + 1.13*) x 1CT 3 , P F a = (-3.04 - 1.27*) x 10" 3 . (31) 



The tree contributions obtained in the three approaches are very close. The signs of the penguin contributions 
are different, and the magnitude in PQCD is larger than those in FA and BBNS: 

|Ppqcd| ~ 5.2|P B bns| , |Ppqcd| ~ 2.7|P FA | . (32) 

The difference in the total penguin contributions is attributed to the treatment of nonfactorizable and annihi- 
lation duagrams. We emphasize that the BBNS results depend on the parameters for factorizable contributions 
(the transition form factors) and the parameters for nonfactorizable contributions (the meson wave funcitons). 
Therefore, it is possible to reverse the sign and to increase the magnitude of the penguin contribution by vary- 
ing these parameters. The penguin contribution in the BBNS approach is smallest because of the cancellation 
between F p and M p . If expressing the amplitudes of the B^ — > tt + tt decay as 

A=V t *\P\e iS -V*\T\, (33) 

with P and T given in Eq. Q) , the strong phases 6 are about 

5 PQCD ~ 123° , 8 BBNS ~ -50° , S FA ~ -113° . (34) 

The predictions for the CP asymmetry of the B° — > 7r ± 7r =F decays are presented in Fig. 5. Because of Eq. (S3), 
the CP asymmetry in PQCD is larger than that (~ 10% at most) in FA and much larger than that (~ 4% at 
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most) in BBNS. Because of Eq. ([34|), the sign of the CP asymmetry in PQCD is opposite to those in FA and 
BBNS. If the B meson wave function in |27|] with lob = 0.4 GeV is adopted, M p and also |Pbbns| will increase 
accordingly, such that the CP asymmetry in the BBNS approach becomes a bit larger. We shall leave the 
resolution of the above discrepancies to future measurements of CP asymmetries. That is, experimental data 
of CP asymmetries will provide an appropriate test of the three approaches. Note that all the three approaches 
give the similar branching ratio B(B^ — ► 7r ± 7r T ) ~ 4.6 x 10~ 6 , which is consistent with the CLEO data p3j. 

As argued in Sec. Ill, to account for the B — > Kir and ttit data simulataneously using the unitarity angle 
4>z = 90°, the ratio of the penguin contribution to the tree contribution must be as large as 0.1. It is easy to 
find 

I^PpqcdI/ITpqcdI = 0.12 , |Fbbns|/|T B bns| = 0.02 , |P F a|/|T fa | = 0.04 , (35) 

from Eq. (]3l|). Due to the smaller ratio in the BBNS and FA approaches, a larger angle 03 > 90° must be 
adopted in order to explain the B — > irn data ||J2^]. As a consequence of the penguin enhancement, the CP 
asymmetry is also large in the PQCD approach. 

VI. SUMMARY 

We summarize the comparision among the FA, BBNS, and PQCD approaches as follows: 

1. The scale independence of predictions in the FA and BBNS approaches is achieved by the generalized 
FA method: a /i-dependent evolution factor is extracted from hadronic matrix elements by considering 
the mixing of the effective operators at a scale below \x. This evolution factor is then combined with 
the /i-dependent Wilson coefficient to form a scale-independent effective Wilson coefficient. However, the 
infrared-cutoff and gauge dependences appear. In the PQCD approach the scale independence is made 
explicit by extending the Wilson evolution down to the hard scale t directly, which is then followed by the 
second-stage RG evolution evaluated based on six-quark amplitudes. Since external quarks are on-shell, 
the gauge invariance of PQCD predictions is guaranteed. The infrared-cutoff dependence, treated as a 
factorization scheme dependence, is removed by the universality of wave functons. 

2. Factorizable contributions (hadronic transition form factors) are not claculable in the FA and BBNS 
approaches, but calculable in the PQCD approach. The inclusion of parton transverse momenta smear 
the end-point singularity, and Sudakov suppression guarantees that almost all contributions come from 
the region with a s /ir < 0.3. 

3. Dynamical penguin enhancement exists in the PQCD approach, but not in the FA and BBNS approaches. 
This enhancement, arising from different evolution effects in the form factors associated with the Wilson 
coefficient a\, and ae, is cruical for the simultaneous explanation of the B — > Kir and 7T7t data using a 
smaller unitarity angle 03 ~ 90° 
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4. Nonfactorizable amplitudes arc neglected in the FA approach, but calculable in the BBNS and PQCD 
approaches. However, without Sudakov suppression, higher-order corrections are not under control in the 
BBNS approach. The end-point singularity needs to be smeared by taking into account parton transverse 
momenta and Sudakov suppression. On the other hand, nonfactorizable amplitudes are real in the BBNS 
approach, but complex in the PQCD approach. The conclusion is the same: nonfactorizable contributions 
are destructive relative to the penguin contributions and negligible for the evaluation of branching ratios. 

5. Annihilation contributions are neglected in the FA and BBNS approaches, but calculable in the PQCD 
appraoch. The helicity suppression applies only to amplitudes associated with O4, but not to those 
associated with 6 ■ It is found that the Oq contribution is of the same order as the penguin contribution, 
and gives a large imaginary part. 

6. Strong phases come from the BSS mechanism and the extraction of the \i dependence from hadronic matrix 
elements in the FA and BBNS approaches, but from annihilation diagrams in the PQCD approach. It has 
been argued that the two sources of strong phases in the FA and BBNS approaches are in fact strongly 
suppressed by the charm mass threshold and by the end-point behavior of meson wave functions. 

7. Because of the large penguin and annihilation contributions in the PQCD approach, the predicted CP 
asymmetry in the B\ — > n ± n T decays dominates over those in the FA and BBNS approaches as shown in 
Fig. 5. Future measurements of CP asymmetries can distinguish the FA, BBNS, and PQCD approaches. 
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Figure Captions 

1. Fig. 1: One of the lowest-order factorizable diagrams for the i?° — > 7r + 7r~ decay. 

2. Fig. 2: Lowest-order nonfactorizable diagrams for the -B° — > 7r + 7r~ decay. 

3. Fig. 3: Lowest-order hard amplitudes for the B® — > tt + it~ decay. 

4. Fig. 4: Diagrams that give strong phases in the FA and BBNS approaches. 

5. Fig. 5: Direct CP asymmetries in the — > 7r ± 7r =F decays. 
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Amplitudes 


Left-handed column 


Right-handed column 


PQCD 


BBNS 


Re(UF T ) 


(a) 


4.39 • 


10~ 2 


(b) 


2.95 • 


10~ 2 


7.34 ■ 10~ 2 


7.57 ■ 10~ 2 


Im(UF T ) 




- 






- 






1.13 • 10~ 3 


Re{f«F p ) 


(a) 


-3.54 


■ 10~ 3 


(b) 


-2.33 


10~ 3 


-5.87- 10~ 3 


-3.04 • 10 -3 


Im(f 7r F p ) 




- 






- 






-1.27 • 10 -3 


Re(f B F a p ) 


(e) 


5.05 • 


10~ 4 


(f) 


-1.94 


10~ 3 


-1.42 • 10~ 3 




Im{f B F p ) 


(e) 


2.19 ■ 


10~ 3 


(f) 


3.72 • 


10~ 3 


5.91 • 10~ 3 




Rp(M T \ 


(c) 


5.02 ■ 


10~ 3 


(d) 


-6.55 


10~ 3 


-1.53 • 10~ 3 


-7.71 • 10~ 4 


Im(M T ) 


(c) 


-3.83 


10~ 3 


(d) 


7.03 • 


10~ 3 


3.20 ■ 10~ 3 




Re(M p ) 


(c) 


-2.29 


10~ 4 


(d) 


2.75 • 


10~ 4 


4.66 ■ 10~ 5 


4.17 ■ 10~ 3 


Im(M p ) 


(c) 


1.95 • 


10~ 4 


(d) 


-3.08 


10" 4 


-1.13 • 10~ 3 




Re(M p ) 


(g) 


1.14 • 


10~ 5 


(h) 


-1.48 


io" 4 


-1.37- 10~ 4 




Im(M p ) 


(g) 


-9.12 


10" 6 


(h) 


-1.27 


io- 4 


-1.36 • 10~ 4 





TABLE I. Amplitudes for the — > ty + ty decay from Fig. 3, where F (M) denotes factorizable (nonfactorizable) 



contributions, P (T) denotes the penguin (tree) contributions, and a denotes the annihilation contributions. Here we 
adopted (f>3 = 90°, Rt = 0.38, and a s (mb) = 0.2552 in the numerical analysis for the BBNS approach. 
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